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Abstract. How much dust can be produced in the early Universe? Does dust 
production depend on the average heavy-metal content of the hosting galaxy? 
Considering supernova explosions, massive stars (Wolf-Rayet, LBV and RSG), 
and relatively massive AGB stars among possible dust-generating objects in the 
early Universe, we find that SN remnants can be regarded as the main source of 
the primordial dust. However, this conclusion is based on highly uncertain (and 
probably over-estimated) dust production rates. Despite all the uncertainties, 
interstellar extinction must be taken into account while observing high-redshift 
objects. 



1. Introduction 

Contemplating the general importance of dust, one should take into account 
the following: 50% of optical radiation emitted since the Big Bang by all 
astrophysical sources has been 'reprocess ed' by dust. Dust is regarded as an 
efficient ISM coolant (Lchncr, Wakkcr &: Savage thus controlling accre- 

tion and profoundly influencing star-formation rates, especially for massive stars 
(e.g. IWolfire fc Cassinelli lll98% . Assembly of H2 on dust grains is far more effi- 
cient than in the gas phase. On cosmological scales, dust may distort the cosmic 
micro wave background and change the far-IR background ( Elfgren Deserti 



I200I). There a re first indic ations of the presence of dust in high-redshift L yman- 
break galaxies ( Ando. et al. 2004 : Ouchi. et al. ,2004, 1 and quasars (Maiolino et al. I 



l2004fL which may have serious im plication for the estimated star-formation rates 
(e.g. ISchiminovich. et al. Il2005t ). 



2. Dust Formation in the Milky Way Galaxy and in the High-Redshift 
Universe 



The detailed galactic census of dust-producing stars ( Gehrz lfl989l l proves that 



in the modern Universe, relatively low-mass stars dominate the scene. All cat- 
egories of evolved, post-main-sequence objects (Miras, carbon stars, late-type 
supergiants, etc.) account for ~ 90% of the stellar dust output, while planetary 
nebulae produce < 1%. Their high-mass counterparts, supernovae and Wolf- 
Rayet stars, amount to < 10%, combined. Dust grains are also manufactured 
via relatively slow-paced accretion in molecular clouds: (1 — 5)x the stellar out- 
put. Apparently, in the 1-2 Gyr-old Universe only the high- mass, M ^ SMq, 
stellar population may be responsible for accumulation of copious amounts of 
dust (Table 1). In the Table 1 we provide: (a) the main categories of dust- 
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producing stars, (b) their initial masses, (c) the total dust yield, per star, (d) 
the chemistry and (e) the size distribution of dust grains, where 'st.' (st andard) 
corresponds to the truncated power-law distribution (M athis et al. lll977.1 . and 
'sm.' (small) is used to emphasize the presence of particles with a -C 1^ sizes. 
We separate the cases of normal, solar metallicity and low-Z environments. 



Table 1. Potential dust-generating stars in the z > 3 Universe 



Category 






Dust 




Ref. 




Mq 


mass/* (Mq) 


composition 


size 




SNt,.^(PopIII) 


140-260 


? ,Z < Zq 


7 


7 




SN II 


> 8 


0.1-0.3, Z -C Zq 

<1, Z - Zq 


Si,C,Fe,Al,Mg,0 
Si,C,Fe,Al,Mg,0 


'St.' (?) 
'st.', ~ 1^ 


1 

2,3 


LBV 


> 75 - 85 

> 30(?) 


? ,Z < Zq 

0.01-0.25, Z ~ Zq 


7 

Si,C,PAH,Al(?) 


7 

~ l/i-|-'sm.' 


4,5,6,7 


WCd 


> 60 - 70 

> 30 


? ,Z < Zq 
10-3 - 10-2, Z ~ Zq 


7 

C(amorph. ) 


7 

~ l/i-|-'sm.' 


8,9,10 


sgB[e] 


7 

> 30 - 60 


? ,Z < Zq 
? ,Z Zq 


7 

Si 


7 

~ iM 


11 


B[e]WD 


> 5 


7 


? 


7 




RSG 


> 8 - 50 

> 8 - 25 


? ,Z < Zq 

10-4 - 10-3, Z ~ Zq 


? 

Si 


? 

~ 0.5/i-l-'st.' 


12,13,14 


AGB 
(OH/IR) 


> 2 - (5 - 6) 

> 2 - 6 


cx Mi„i{Z/ZQ),Z < Zq 
SJ 10~^Mi„i,Z ~ Zq 


7 

Si, C, ice 


7 

'St.'' (?) 






iMaiolino et al. I (E)ni). 12]^. 

IWaters et a]_ I (Il997l). 



Clayton et ah 1 (Exvl). f3l ISugerman et a~ 
[51 IVoorsetal. I (Il999l). [fil lyoors et al 



[41 IWaters et a1 I (Il997l). f5l IVoors et al. I (Il999l). [fil 
[71 iKing et ah I (120021). [81 IWillia.ms. van der Hucht fcThe i iim?SYii. iy> 
al. I (1199^, [10 1 iMarchenko et al I (|2002l), [ lll.Molster et al. I (|2002l ) 



iMathis et _ , 

[12] ISeab fc Snow I (Il989l ). [13] I.Tural (|l99fil ). [14] .Smith et al. I (1200 j ) 



Judging by the multitude of question marks, not much is known about 
dust formatio n in low-Z environments. The pair-instability supernovae (SN^-y: 
iHeger fc Wo~slev , .20o3 ) belong to the broad category of Population III ob- 
jects. Though their general characteristics and ability to produce dust are yet 
to be established, one may assume their dust yields to be comparable with 
the present-epoch SN events, thus making them the major dust sources in the 
high-redshift universe. However, one should notice the substantial difference 
between the total dust outputs of SN events provi ded bv different rese arch 
groups: the re latively high, M^ust ~ l-^O; estimates of lDunne et al. I (120031) and 
Morgan et al. (2003), in line with theoretic al expectations (iTodini &: Ferrara. I 



)rga 

ml 
IE 



200lfl. vs. th7Mw , ,o^ ^ IQ-^Mg, values from lDwek et al. I (|l992l ). iDwek I (|2004l ) 



and lPozzo et al. I (|2004l ). There is a strong indicatio n that at least in so me SN 
event s (the Crab nebula, SN2002hh and SN2002ic: Gre en et al. 2004; Barlo w et ai~l 
l200,l iKotak et al. 11200.^ dust comes from a progenitor, either a luminous blue 
variable (LBV), red supergiant (RSG) or carbon-rich Wolf-Rayet star (WCd). 
With the lower yields for SN events, the integral output from RSG and WCd 
stars may rival the SN categ ory. The relatively rare B supergiants with for- 
bidden emission lines, sgB[e] ( Lamers et al. iri998r ). as well as newly discovered 
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(and probably associated with sgB je]) class of early-type, lumi nous stars with 
warm circumstellar dust, B[e]WD ((Miroshnichenko et al. f2005 l). could be con- 
sidered as minor contributors, unless the number of related sources is grossly 
under-estimated. One more category emerges at z < 9 and gradually comes to 
a complete dominance over the dust production in the modern universe: the 
asymptotic giant-bran ch stars (AGB). This group shows a clear dependence of 
dust production on Z (Ivan Loon II2000D . This dependence is less pronounced for 
SNe (iTodini fc Ferrara 1120011 ^" For the remaining categories of massive dust- 
producing stars one may assume that dust produc tion depends on the mass 
loss ra te and use the general M=f(Z) relationship of IVink. de Koter Sz Lamers I 

(l2nnih . 



3. Interstellar Extinction for High-Redshift Objects 

Estimating an average e xtinction for high- r edshif t objects, we adopt and slightly 
modify the approach of lLoeb Sz Haimanl (^O^)- Namely, the dust absorption 
coefficient, a^, is expressed as 

au{z, Z) = pdustiz, Z)k^{Z), 

with dust opacity k,v{Z) represented by the Galactic law ((Mathis 'l99[f) for a_Z~ 
Z(7^ e nvironment, or the Small Magellanic Cloud dependence (Cartledge e t al. I 
l2005l ^ for Z < Zq; z defines the redshift and Z denotes the metallicity. The dust 
density takes the form 

Pdustiz, Z) = QbPc{l + zf^Fi{z)fdep,i{z)fdust,ii^), 

i 

where the sum runs over different categories of dust producers (Table 1), Pc = 
9.7 g cm~^ provides the current critical density of the universe, and fif, = 
0.044 gives the total baryonic density. The mass fraction of dust, fdust,i{Z),, 
deposited by a given star depends on the category of dust-producing stars and 
the ambient metallicity, while the mass fraction of stars being able to produce 
dust is calculated as 

M2i Mu{z) 
MU Ma{z) 

with Md{z), Mu{z) and x depending on the redshift, and Mlj, M2j provided by 
Table 1. Hence, the variable x could be anywhere between x=0.5 ('top-heavy' 
mass function) and 1.35 (classical form) at z > 10, then converging to x=1.35 
for z < 10. We define the Fi(z) term as 

20 

Fiiz) = j m{z)-^exp{ — ^-^^)fstar{z')dz' , 



where the efficiency of star formation fstar{z) is lowered by the presence of com- 
pact objects (neutron stars, black holes, white dwarfs) in the overall 'recycling' 
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loop: fL.Jz) = f.,nr(z)(l - eL,nJz + dz)) , with ^ 
taken froml Prorv et alTI (|2005l ). The upper h mit, z= 



= 0.10-0.15 and fstar{z) 
20, is imposed by the evo- 
lution of Population III objects (e.g. Chou dhurv fc: Ferrara . .2005 ) . We adopt 
the mass fraction of bar yons assembled into collapsed objects, Fcoi{z), following 
iHaiman Sz Loeb 1 (^Q^). The exp{—{tz — tz')/T) te rm provides the dust sur- 
vival timescales, with T ranging from 0.1 to 1.0 Gyr (iDraine fc Salneter 1979) . 
The term r]i{z) introduces evolutionary timescales (i.e. cutoffs) for different cat- 
egories of dust producers: rii{z) = 1 for z < Zcr,i, and 7?i(z) = for z > Zcr,i- 
Then, the optical depth of a dusty medium at a given redshift > 3 is 



Tdust{t^, Zs, Z) 



c 



a 



v{l+z) 



{Z,Z) 



(l+z)5/2 



dz, 



with H n and c as universal constants. Following the arguments of lLoeb k, Haiman I 
(1993), we ignore dust production at z < 3, as it will be generally confined to 
individual galaxies (dominance of AGB stars; hence, rather low velocities of dust 
ejecta) rather than ~homogeneously distributed along the line of sight. 




I I I I 
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Figure 1. The dust opacity for a z=6 object with variable yield from the 
SN ejecta. 



Exploring the d i fferen t parameters of the model, we confirm the conclusion 
of iLoeb fc Haiman I (|l997l ) that dust chemistry has relatively small influence 
on the overall opacity. More serious is the unce rtainty in the dust yield s of 
SNe. Considering the theoretical predictions fe.g. iTodini &: Ferrara ll200lh . we 
find that they result in a,n incLppropriatcly high. T^ust 

for a z=6 object (Fig. 1). 

Lowering them by an order of magnitude brin gs the theoretical yie lds closer to 
some estimates from recent observations (e.g. iBarlow et al. II2005I ). With the 
appropriately adjusted SN output, the average dust opacity may be neglected 
for 2: ~ 5 objects, unless they pass through a star-burst episode, and should be 
taken into account for z ^ 5 objects (Fig. 2). 

Grouping the dust-producing stars into 3 general categories: SN events vs. 
massive stars (RGB, LBV, WCd) vs. AGB, and varying the yield from the 
SN stars, as well as adjusting the dependence of the dust production rate on 
Z (metallicity) for the massive stars, we plot two extreme scenarios in Fig. 3. 



5 




Considering the shares of dust production in the modern universe, one may ex- 
pect that the model with the lowest SN yield and absence of a steep dependence 
of the dust production rate in massive stars on the ambient metallicity (right 
panel of Fig. 3) may be closer to reality. 

The calculations also show that, on average, the 'survival' timescales of the 
primordial dust, T < 410^ yr, closely match the current-epoch expectations for 
the Galaxy (I.Tones et al. 1119941) . 




10 a 6 4 10 8 6 4 



z z 

Figure 3. Three main groups of dust-producing stars: SN (full lines), mas- 
sive stars (RSG-hLBV-l-WCd: dotted lines) and AGB (dashed lines). Left 
panel: the SN rate is lowered to 1/10 of the theoretical predictions and the 
dust production rate for massive stars depends on Z. Right panel: SN rate 
=1/30 theoretical, no dependence on Z for massive stars. 



4. Conclusions 

• SN events should be regarded as a major source of dust in the high-redshift 
(z > 3) universe. This conclusion may be independent of the source of dust, 
unless different sources provide substantially different extinction curves 
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(e.g. iMaiolino et al. Il2004^ . It could be either primordial dust coming 
from a progenitor (LBV, RSG or WCd star), or dust produced in the SN 
ejecta. However, the most important issue is the ability of dust to survive 
in the hostile environment: shocks, UV radiation. On an optimistic note, 
one may mention the case of WCd stars where, facing a similar challenge as 
in the shocked environments of SN ejecta, t he grains of amorphous carbon 
manage to survive for at least 10 years (|Marchenko et al. ll2n?)l ). thus 
effectively reaching (and enriching) the ISM. 

• Our calculations show that, in order to b e comparable to the known output 
of SN events in the Galaxy ((Gehrz Ill989l ^ , the theoretical estimates of dust 
production in SNe populating the 1-2 Gyr-old universe should be lowered 
by an order of magn itude, thus providing a better match to the recent (e.g. 
iBarlow et al. 11200.^ observations. 

• Overall, IS extinction should be appropriately taken into account for all 
z ^ 5 objects in order to estimate their true properties, especially realizing 
that for some of them, due to the enhanced star-formation rate, the local 
extinction may substantially exceed the average values (cf. Fig. 2). 
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